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ABSTRACT 

The  astrometric  calibration  of  the  Sloan  Digital  Sky  Survey  is  described.  For  point  sources  brighter  than 
r  ~  20,  the  astrometric  accuracy  is  45  mas  rms  per  coordinate  when  reduced  against  the  USNO  CCD  Astro- 
graph  Catalog  and  75  mas  rms  when  reduced  against  Tycho-2,  with  an  additional  20-30  mas  systematic  error 
in  both  cases.  The  rms  errors  are  dominated  by  anomalous  refraction  and  random  errors  in  the  primary  refer¬ 
ence  catalogs.  The  relative  astrometric  accuracy  between  the  r  filter  and  each  of  the  other  filters  (w,  g ,  z,  z)  is 
25-35  mas  rms.  At  the  survey  limit  (r  ~  22),  the  astrometric  accuracy  is  limited  by  photon  statistics  to 
approximately  100  mas  rms  for  typical  seeing.  Anomalous  refraction  is  shown  to  contain  components 
correlated  over  2°  or  more  on  the  sky. 
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1.  INTRODUCTION 

The  Sloan  Digital  Sky  Survey  (SDSS)  is  obtaining  images 
of  one-quarter  of  the  sky  in  five  broad  optical  bands  (zz,  g,  r, 
z,  and  z;  Fukugita  et  al.  1996;  Smith  et  al.  2002;  Hogg  et  al. 
2001),  with  95%  completeness  limits  for  point  sources  of 
22.0,  22.2,  22.2,  21.3,  and  20.5  mag,  respectively.  The  survey 
will  produce  a  database  of  roughly  108  galaxies,  108  stars, 
and  106  quasars  with  accurate  photometry,  astrometry,  and 
object  classification  parameters  (York  et  al.  2000).  Spectros¬ 
copy,  covering  the  wavelength  range  of  3800-9200  A  at 
R  «  1800,  will  be  obtained  for  roughly  1,000,000  galaxies, 
100,000  quasars,  and  another  50,000  stars  and  serendipitous 
objects.  (Currently,  funding  is  available  to  complete 
approximately  two-thirds  of  the  survey  area.)  At  the  time  of 
this  writing,  approximately  3600  deg2  of  sky  have  been 
imaged,  and  spectra  have  been  obtained  for  400,000  objects. 

This  paper  describes  the  astrometric  calibration  of  the 
SDSS.  Accurate  placement  of  fibers  for  the  survey  spectro¬ 
scopy  requires  an  accuracy  for  astrometry  of  180  mas  rms 
per  coordinate  (all  errors  and  accuracies  quoted  in  this 
paper  are  rms  per  coordinate).  In  fact,  the  final  accuracy  is 
considerably  better.  Depending  on  the  reference  catalog 
available  in  a  given  area  of  sky,  for  point  sources  brighter 
than  r  ~  20,  the  astrometry  is  accurate  to  45  mas  rms  when 
reduced  against  the  USNO  CCD  Astrograph  Catalog 
(UCAC;  Zacharias  et  al.  2000)  and  75  mas  rms  when 
reduced  against  Tycho-2  (Hog  et  al.  2000),  with  additional 
systematic  errors  of  20-30  mas  in  both  cases.  The  accuracy 
of  the  relative  astrometry  between  the  r  detection  and  the  u, 
g,  z,  and  z  detections  is  typically  25-35  mas. 


Section  2  describes  the  layout  of  the  imaging  camera, 
which  motivates  much  of  the  astrometric  calibration  pro¬ 
cess.  Section  3  gives  an  overview  of  the  imaging  data  proc¬ 
essing  pipelines,  providing  a  context  for  the  detailed 
discussion  of  the  astrometric  pipeline.  Section  4  discusses 
centroiding  algorithms.  The  astrometric  pipeline  and  cali¬ 
bration  procedures  are  described  in  §  5,  and  the  accuracy  of 
the  resulting  astrometry  is  assessed  in  §  6.  Differences 
between  the  calibrations  described  here  and  those  employed 
in  previously  released  SDSS  data  are  described  in  §  7.  A  few 
concluding  remarks  are  given  in  §  8. 

2.  IMAGING  CAMERA 

The  imaging  camera  (Gunn  et  al.  1998;  Fig.  1)  consists  of 
54  CCDs  in  eight  Dewars  and  spans  2? 3  on  the  sky.  Thirty 
of  these  CCDs  are  the  main  imaging/photometric  devices, 
each  a  SITe  (Scientific  Imaging  Technologies,  formerly  Tek¬ 
tronix)  device  with  20482  24  jum  square  pixels.  They  are 
arranged  in  six  column  Dewars  of  five  CCDs  each,  one 
CCD  for  each  filter  bandpass  (u,  g,  r,  z,  z)  in  each  column. 
The  camera  is  operated  in  TDI  (time  delay  and  integrate), 
or  scanning,  mode,  for  which  the  telescope  is  driven  at  a  rate 
synchronous  with  the  charge  transfer  rate  of  the  CCDs. 
Objects  on  the  sky  drift  down  columns  of  the  CCDs  so  that 
five-color  photometry  is  obtained  over  an  interval  of  5 
minutes  with  no  dead  time  for  shuttering  and  readout.  The 
effective  integration  time  is  53.9  s  at  the  chosen  (sidereal) 
scanning  rate  (the  pixel  scale  is  0'/396  pixel-1). 

The  other  24  CCDs  in  two  additional  Dewars  are  also 
SITe  chips  of  width  2048  24  jum  pixels,  but  they  have  only 
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Fig.  1. — Layout  of  the  imaging  camera.  The  large  colored  squares  represent  the  photometric  CCDs,  five  in  each  of  six  Dewars,  which  run  vertically  on  the 
page.  The  smaller  rectangles  above  and  below  show  the  22  astrometric  and  two  focus  CCDs  in  two  Dewars  running  horizontally.  The  arrows  indicate  the 
direction  of  motion  of  stars  across  the  camera  during  a  scan. 


400  rows  in  the  scanning  direction.  These  Dewars  are  ori¬ 
ented  perpendicular  to  the  photometric  Dewars,  and  one 
lies  across  the  top  of  the  columns,  the  other  across  the  bot¬ 
tom.  Each  Dewar  contains  12  CCDs.  Two  of  these  CCDs 
(one  in  each  Dewar)  are  used  to  dynamically  determine 
changes  in  focus.  The  other  22  CCDs  are  used  to  tie  obser¬ 
vations  to  an  astrometric  reference  frame.  They  saturate  at 
brighter  magnitudes  than  the  photometric  CCDs  because 
they  (1)  have  a  shorter  integration  time  (10.5  s),  since  they 
have  400  versus  2048  rows  in  the  scanning  direction,  (2)  are 
behind  a  neutral  density  filter  of  3.0  stellar  magnitudes,  and 
(3)  have  somewhat  poorer  quantum  efficiency.  The  astrome¬ 
tric  CCD  filter  bandpass  is  similar  to  the  photometric  r 
bandpass,  although  it  extends  ~400  A  redward  of  the  r  red 
cutoff.  The  astrometric  CCDs  saturate  at  r  ~  8  and  thus  can 
utilize  all  but  the  brightest  10%  of  the  astrometric  standard 
stars  in  Tycho-2.  The  photometric  r  CCDs  saturate  at 
r  ~  14,  while  good  centroids  can  be  measured  on  the  astro- 
metric  CCDs  for  stars  as  faint  as  r  ~  17.  Hence,  there  is  an 
overlap  of  about  3  mag  between  faint  stars  with  good  cent¬ 


roids  on  the  astrometric  chips  and  bright  unsaturated  stars 
on  the  r  chips,  enabling  the  transfer  of  catalog  stars  detected 
on  the  astrometric  CCDs  to  the  r  CCDs  and  thus  the  astro- 
metric  calibration  of  stars  detected  on  the  r  CCDs  against 
bright  reference  catalogs.  Similarly,  there  is  an  overlap  of  at 
least  5  mag  between  stars  with  good  centroids  detected  on 
the  r  CCDs  and  the  other  photometric  CCDs,  allowing  the 
r  astrometric  calibrations  to  be  transferred  to  the  other 
photometric  CCDs. 

Twelve  of  the  astrometric  CCDs  are  aligned  with  the 
six  columns  of  photometric  CCDs,  with  an  astrometric 
chip  leading  and  trailing  each  column.  The  other  10 
astrometric  chips  bridge  the  gaps  between  the  columns. 
Hence,  the  astrometric  chips  span  the  full  2?3  width  of 
the  camera.  As  outlined  below,  the  bridge  astrometric 
CCDs  are  not  used  in  the  current  version  of  the  astrome¬ 
tric  pipeline.  They  were  originally  intended  to  account 
for  differential  shifts  between  the  photometric  Dewars, 
but  the  camera  has  proven  to  be  extremely  stable,  and 
the  technique  was  not  needed. 


No.  3,  2003 


SDSS  ASTROMETRIC  CALIBRATION 


1561 


3.  DATA  PROCESSING  OVERVIEW 

The  large  volume  of  data  generated  by  the  SDSS  requires 
a  set  of  automated  data  processing  pipelines  to  acquire  and 
archive  the  raw  data,  detect  and  measure  object  parameters, 
astrometrically  and  photometrically  calibrate  the  objects, 
select  spectroscopic  targets,  design  the  spectroscopic  plates, 
extract  spectra,  and  measure  various  spectral  parameters. 
The  end-to-end  processing  of  SDSS  data  is  discussed  in 
detail  in  Stoughton  et  al.  (2002).  Those  parts  of  data  pro¬ 
cessing  that  are  relevant  to  the  astrometric  calibrations  are 
briefly  summarized  here  to  give  context  to  the  detailed  dis¬ 
cussion  of  the  astrometric  calibrations. 

Each  imaging  drift  scan  is  processed  serially  through  the 
following  pipelines: 

1 .  The  real-time  data  acquisition  system  (DA;  Petravick 
et  al.  1994)  acquires  the  raw  data  as  they  are  read  off  the 
CCDs,  buffers  them  to  disk,  archives  them  to  tape,  and  per¬ 
forms  real-time  quality  analysis  (it  is  the  only  real-time  com¬ 
ponent  of  data  processing).  The  data  stream  from  each 
CCD  is  broken  up  into  frames  (2048  columns  by  1361  rows) 
to  facilitate  data  processing.  DA  finds  bright  stars  on  the 
frames  from  the  astrometric  CCDs,  cuts  out  a  29  x  29  pixel 
(1 1  '/5  x  11 '.'5)  subraster  (“postage  stamp”)  centered  on 
each  star,  and  measures  parameters,  including  a  centroid 
and  flux,  for  each  star  (the  astrometric  frames  are  not 
saved). 

2.  The  serial  stamp  collection  pipeline  (SSC;  Lupton  et 
al.  2001)  cuts  out  postage  stamps  and  measures  centroids 
for  bright  stars  on  the  photometric  CCD  frames.  It  mea¬ 
sures  stars  at  the  expected  positions  based  on  the  DA  detec¬ 
tions  on  the  astrometric  CCDs  in  the  same  camera  column 
and  includes  a  star  finder  to  detect  additional  stars  not 
detected  on  the  astrometric  CCDs. 

3.  The  postage  stamp  pipeline  (PSP;  Lupton  et  al.  2001) 
uses  postage  stamps  for  bright  stars  cut  out  by  SSC  to  meas¬ 
ure  the  point-spread  function  (PSF)  as  a  function  of  posi¬ 
tion  across  each  frame  (for  the  photometric  CCDs  only). 
The  star  centroids  for  all  postage  stamps  are  remeasured 
using  an  algorithm  that  corrects  for  asymmetries  in  the 
PSF,  as  described  in  §  4. 

4.  The  astrometric  pipeline  (Astrom)  matches  up  the 
stars  detected  by  DA  for  the  astrometric  CCDs,  and/or 
those  detected  by  SSC  for  the  photometric  CCDs  with  cor¬ 
rected  centroids  produced  by  PSP,  to  a  primary  astrometric 
reference  catalog.  Using  these  matches,  it  derives  a  separate 
set  of  transformation  equations  for  each  frame,  converting 
frame  row  and  column  to  J2000.0  catalog  mean  place  celes¬ 
tial  coordinates.  This  procedure  provides  the  astrometric 
calibration  for  the  SDSS  and  is  the  topic  of  this  paper.  The 
pipeline  itself  is  described  in  detail  in  §  5. 

5.  The  imaging  pipeline  (Frames;  Lupton  et  al.  2001) 
corrects  the  raw  frames  for  instrumental  effects,  detects 
objects  on  the  corrected  frames,  matches  up  detections  of 
the  same  object  in  different  filters,  deblends  overlapping 
objects,  and  measures  parameters  for  the  objects.  It  uses 
the  astrometric  calibrations  produced  by  Astrom  to  pre¬ 
cisely  match  up  the  separate  detections  in  the  different  fil¬ 
ters  for  each  object.  This  requires  that  Astrom  run  before 
Frames  and  explains  why  the  astrometric  calibrations  are 
based  on  centroids  measured  by  PSP  and  DA,  rather 
than  on  Frames’s  centroids.  Frames  and  PSP  use  the 
same  centroiding  algorithm.  The  final  SDSS  object  posi¬ 
tions  are  derived  by  applying  the  astrometric  calibrations 


produced  by  Astrom  to  the  centroids  measured  by 
Frames. 

In  summary,  the  astrometric  calibrations  are  produced  by 
Astrom,  processing  a  single  drift  scan  at  a  time,  using  cen¬ 
troids  measured  by  DA  and/or  PSP  and  producing  a  set  of 
transformation  equations  for  each  frame  that  are  applied  to 
centroids  measured  by  Frames. 


4.  CENTROIDS 

An  object’s  centroid  is  defined  as  the  first  moment  of  its 
light  distribution.  Since  this  is  a  noisy  estimate  for  most 
objects  in  the  survey,  the  following  technique  is  used  to  bet¬ 
ter  estimate  the  centroid.  First,  an  object’s  image  is 
smoothed  using  a  two-dimensional  Gaussian  with  an  adap¬ 
tive  smoothing  length  scaled  to  the  PSF  in  that  frame. 
Quartic  interpolation  is  used  to  find  the  maximum  in  a  3  x  3 
pixel  subraster  centered  on  the  peak  pixel  in  the  smoothed 
image.  This  gives  a  biased  estimate  of  the  centroid  in  the 
presence  of  an  asymmetric  PSF.  PSP  uses  bright  stars  to 
determine  the  shape  of  the  PSF  as  a  smoothly  varying  func¬ 
tion  of  CCD  column  and  row  for  each  frame.  Thus,  the  PSF 
at  the  position  of  each  object  is  determined  to  high  accuracy. 
The  centroid  of  the  PSF  at  the  position  of  the  object  is  mea¬ 
sured  using  both  the  first  moment  and  quartic  interpolation 
after  smoothing  by  the  PSF.  The  difference  in  centroids 
measured  for  the  PSF  using  the  two  algorithms  is  a  measure 
of  the  bias  introduced  by  the  use  of  the  quartic  interpolation 
algorithm.  This  bias  is  typically  a  few  tens  of  milliarcsec- 
onds  but  can  be  as  large  as  100  mas.  The  difference  is  added 
to  the  quartic  interpolation  centroid  measured  for  the 
object,  yielding  a  high  signal-to-noise  ratio  estimate  of  the 
first-moment  centroid  of  the  object.  A  more  detailed  discus¬ 
sion  of  the  centroiding  algorithm  will  appear  in  the  forth¬ 
coming  paper  on  SDSS  photometry  (R.  Lupton  et  al.  2003, 
in  preparation). 

PSP  and  Frames  both  use  the  bias-corrected  centroiding 
algorithm  described  above.  This  is  necessary  because  the 
astrometric  calibrations  are  produced  from  centroids  mea¬ 
sured  by  PSP  but  are  applied  to  centroids  measured  by 
Frames.  Since  neither  DA  nor  SSC  has  an  estimate  of  the 
PSF,  they  measure  uncorrected  quartic  interpolation  cent¬ 
roids,  using  a  constant  smoothing  length  of  cr  =  1.2  pixels 
(~0'/48).  When  reducing  against  Tycho-2,  Astrom  uses 
DA’s  uncorrected  centroids  for  objects  detected  on  the 
astrometric  CCDs;  however,  the  biases  thus  introduced  are 
largely  removed  in  the  calibrations  and  are  small  compared 
to  more  dominant  atmospheric  effects. 

In  addition  to  its  row  and  column  centroid,  each  detec¬ 
tion  of  a  star  is  characterized  by  the  terrestrial  time  (TT) 
when  the  star  was  at  midexposure,  its  flux  on  that  frame, 
and  an  instrumental  color  based  on  detections  on  two  CCDs 
in  the  same  camera  column  during  the  same  drift  scan  (r—i 
for  the  r,  z,  z,  and  astrometric  CCDs,  u—g  for  the  u  CCDs, 
and  g—r  for  the  g  CCDs).  Astrom  converts  instrumental 
fluxes  and  colors  to  approximate  calibrated  magnitudes  and 
colors  using  a  separate  zero  point  for  each  CCD  (derived 
from  a  recent  photometric  calibration),  the  air  mass  at  the 
middle  of  the  frame,  and  average  values  for  the  extinction 
(0.520  mag  per  unit  air  mass  in  u,  0.200  in  g ,  0.120  in  r, 
0.080  in  z,  and  0.065  in  z).  These  calibrated  magnitudes  are 
used  by  Astrom  for  quality  analysis  and  to  reject  faint  stars. 
The  calibrated  colors  are  used  for  differential  chromatic 
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refraction  (DCR)  calculations.  The  expected  errors  in  the 
approximate  magnitudes  and  colors  do  not  contribute  sig¬ 
nificantly  to  the  overall  astrometric  error  budget.  Only 
unsaturated  stars  brighter  than  20th  magnitude  are  used  to 
derive  the  astrometric  calibrations. 

5.  CALIBRATIONS 

The  r  CCDs  serve  as  the  astrometric  reference  CCDs  for 
the  SDSS.  That  is,  the  positions  for  SDSS  objects  are  based 
on  the  r  centroids  and  calibrations.  One  of  two  reduction 
strategies  is  employed,  depending  on  the  coverage  of  astro- 
metric  catalogs: 

1 .  Whenever  possible,  stars  detected  on  the  r  photometric 
CCDs  are  matched  directly  with  stars  in  UCAC.  UCAC 
extends  down  to  R  =  16,  giving  approximately  2-3  mag  of 
overlap  with  unsaturated  stars  on  the  photometric  r  CCDs. 
The  astrometric  CCDs  are  not  used. 

2.  If  UCAC  coverage  is  not  available  to  reduce  a  given 
imaging  scan,  detections  of  Tycho-2  stars  on  the  astrometric 
CCDs  are  mapped  onto  the  r  CCDs  (all  Tycho-2  stars  satu¬ 
rate  on  the  r  CCDs)  using  bright  stars  that  have  sufficient 
signal-to-noise  ratio  on  the  astrometric  CCDs  and  are  unsat¬ 
urated  on  the  r  CCDs.  (This  mode  was  originally  expected 
to  be  the  norm;  at  the  time  the  camera  was  designed,  UCAC 
was  several  years  away  from  beginning  observations.) 

The  r  CCDs  are  therefore  calibrated  directly  against  the 
primary  astrometric  reference  catalog.  Frames  uses  the 
astrometric  calibrations  to  match  up  detections  of  the  same 
object  observed  in  the  other  four  filters.  The  accuracy  of  the 
relative  astrometry  between  filters  can  thus  significantly 
impact  Frames,  in  particular  the  deblending  of  overlapping 
objects,  photometry  based  on  the  same  aperture  in  different 
filters,  and  detection  of  moving  objects.  To  minimize  the 
errors  in  the  relative  astrometry  between  filters,  the  u ,  g ,  z, 
and  z  CCDs  are  calibrated  against  the  r  CCDs. 

Each  drift  scan  is  processed  separately.  All  six  camera  col¬ 
umns  are  processed  in  a  single  reduction.  In  brief,  stars 
detected  on  the  r  CCDs  if  calibrating  against  UCAC,  or 
stars  detected  on  the  astrometric  CCDs  transformed  to  r 
coordinates  if  calibrating  against  Tycho-2,  are  matched  to 
catalog  stars.  Transformations  from  r  pixel  coordinates  to 
catalog  mean  place  (CMP)  celestial  coordinates  are  derived 
using  a  running-mean  least-squares  fit  to  a  focal  plane 
model,  using  all  six  r  CCDs  together  to  solve  for  both  the 
telescope  tracking  and  the  r  CCDs’  focal  plane  offsets,  rota¬ 
tions,  and  scales,  combined  with  smoothing  spline  fits  to  the 
intermediate  residuals.  These  transformations,  comprising 
the  calibrations  for  the  r  CCDs,  are  then  applied  to  the  stars 
detected  on  the  r  CCDs,  converting  them  to  CMP  coordi¬ 
nates  and  creating  a  catalog  of  secondary  astrometric  stan¬ 
dards.  Stars  detected  on  the  u,  g,  z,  and  z  CCDs  are  then 
matched  to  this  secondary  catalog,  and  a  similar  fitting  pro¬ 
cedure  (each  CCD  is  fitted  separately)  is  used  to  derive 
transformations  from  the  pixel  coordinates  for  the  other 
photometric  CCDs  to  CMP  celestial  coordinates,  compris¬ 
ing  the  calibrations  for  the  u,  g,  z,  and  z  CCDs.  The  remain¬ 
der  of  this  section  describes  the  calibration  procedure  in 
detail. 


5.1.  Great  Circle  Coordinates 

In  order  to  minimize  curvature  and  time  differences  as 
objects  drift  across  the  focal  plane,  all  SDSS  drift  scans  are 


conducted  along  CMP  J2000.0  great  circles.  Each  such  great 
circle  (referred  to  as  a  44  stripe  ”)  is  defined  by  the  inclination 
and  the  right  ascension  of  the  ascending  node  where  the 
great  circle  crosses  the  J2000.0  celestial  equator  (Stoughton 
et  al.  2002  provides  more  information  on  the  various  coordi¬ 
nate  systems  employed  in  the  survey).  The  “boresight”  is 
defined  as  the  position  in  the  focal  plane  that  nominally 
tracks  the  great  circle.  A  minimum  of  two  scans  is  required 
along  each  great  circle,  one  with  the  boresight  offset  22.74 
mm  below  (perpendicular  to  the  scan  direction)  the  camera 
center  (referred  to  as  the  north  4  4  strip  ”  of  the  stripe)  and 
one  with  the  boresight  offset  22.74  mm  above  the  camera 
center  (referred  to  as  the  south  strip).  The  two  strips  are 
interlaced  to  fill  in  the  gaps  between  the  camera  columns, 
such  that  two  strips  together  cover  a  swath  about  2?53  wide 
centered  along  the  great  circle.  The  offsets  between  strips  are 
chosen  so  that  the  individual  scan  lines  from  each  of  the  six 
CCD  columns  overlap  slightly  (^0.'5)  with  adjacent  scan 
lines  when  interlaced.  Adjacent  stripes,  separated  by  2?5, 
also  overlap  at  their  edges.  The  boresight  tracking  rate  is 
adjusted  by  the  telescope  control  computer  (TCC)  to  be 
constant  in  observed  place. 

Reductions  are  simplified  by  working  in  a  coordinate  sys¬ 
tem  in  which  the  tracked  great  circle  is  the  equator  of  the 
coordinate  system.  In  this  coordinate  system,  referred  to 
throughout  as  the  44  great  circle  coordinate  system,”  the  lati¬ 
tude  of  an  observed  star  never  exceeds  about  1?3;  thus,  the 
small-angle  approximation  can  be  used,  and  lines  of  con¬ 
stant  longitude  are  to  a  high  approximation  perpendicular 
to  lines  of  constant  latitude.  Longitude  and  latitude  in  great 
circle  coordinates  are  referred  to  with  the  symbols  m  and  z/, 
respectively.  Along  the  great  circle,  v  =  0,  m  increases  in  the 
scan  direction,  and  the  origin  of  m  is  chosen  so  that 
j u  =  (^2000  at  the  ascending  node  (where  the  great  circle 
crosses  the  J2000.0  celestial  equator).  The  conversion  from 
great  circle  coordinates  to  J2000.0  celestial  coordinates  is 
then 


tan(a2ooo  -  Mo) 


sin (ju  —  ju0)  cos  v  cos  z  -  sin  v  sin  z 

COS (jLL  ~  JUq)  cos  v 


sin(<52ooo)  =  sin(//  -  Mo)  cos  v  sin  z  +  sin  v cos  z  , 


(1) 

(2) 


where  z  and  Mo  are  the  inclination  and  J2000.0  right  ascen¬ 
sion  of  the  great  circle  ascending  node,  respectively 
(ju0  =  95°  for  all  survey  stripes).1 


5.2.  Pseudo-Catalog  Place 

As  stated  above,  the  telescope  boresight  tracks  (or  should 
track)  a  J2000.0  great  circle  in  CMP.  The  distance  of  a  star’s 
image  from  the  boresight  on  the  telescope  focal  plane  is  pro¬ 
portional  to  the  difference  in  observed  place  between  the 
star  and  the  boresight  pointing.  The  reductions  are  there¬ 
fore  performed  in  pseudo-catalog  place  (PCP),  a  non¬ 
standard  place  defined  as 

^*,pcp  =  Mb, cmp  +  (A*, op  —  Mb, op)  5  (3) 

^*,PCP  =  Vb, CMP  +  (^*,OP  -  Vb, op)  ?  (4) 

where  the  subscripts  *  and  B  indicate  the  star  and  boresight, 
respectively,  and  the  subscripts  PCP,  CMP,  and  OP  signify 


1  A  few  early  commissioning  scans  and  some  ongoing  calibration  scans 
track  great  circles  that  are  not  survey  stripes. 
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pseudo-catalog  place,  catalog  mean  place,  and  observed 
place,  respectively  (for  definitions  of  various  astrometric 
places,  see  Hohenkerk  et  al.  1992).  Local  place  could  just  as 
well  have  been  used;  however,  using  PCP  facilitates  the  gen¬ 
eration  of  quality  analysis  information  regarding  the  tele¬ 
scope  tracking  and  camera  rotation  over  the  course  of  a 
scan.  That  is,  the  pipeline  explicitly  solves  for  parameters 
describing  the  actual  telescope  scanning. 

5.3.  Calibration  Equations 

While  the  reductions  are  performed  in  PCP,  ultimately  a 
transformation  from  pixel  coordinates  to  CMP  is  required. 
Throughout  all  SDSS  data  processing,  the  data  from  each 
CCD  in  a  drift  scan  are  broken  up  into  contiguous  frames  of 
2048  columns  (columns  are  parallel  to  the  scan  direction)  by 
1361  rows  (perpendicular  to  the  scan  direction).  Astro- 
metric  calibrations  are  generated  as  a  separate  set  of 
equations  for  each  frame,  converting  frame  row  (x),  frame 
column  (y),  and  star  color  to  CMP  great  circle  coordinates 


(/Amp?  *Amp): 
for  r  —  i  <  (r  —  z)0  : 

x'  =  X  +  go  +  giy  +  g2y 2  +  giy2  +  px color  ,  (5) 

/  =  y  +  K  +  hy  +  h2y2  +  h3y2  +  py color  ;  (6) 

for  r  —  i  >  (r  —  z)0  : 

x'  =  X  +  go  +  giy  +  g2y2  +  giy2  +  qx  ,  (7) 

/  =  y  +  ho  +  hiy  +  h2y2  +  hy2  +  qy  ;  (8) 

y«CMP  =  a  +  bx'  +  cy' ;  (9) 

z'cmp  =  d  +  ex'  +fy'  .  (10) 


The  transformation  from  (x,  y)  to  (x',  /)  corrects  for  optical 
distortions  (which,  in  TDI  mode,  are  a  function  of  column 
only)  and  DCR.  For  u  and  g  frames,  DCR  is  modeled  as  a 
linear  function  of  color  (u—g  for  u  frames,  g—r  for  g  frames) 
for  stars  bluer  than  (r— z')0  =  1.5,  and  a  constant  for  redder 
stars.  For  r,  z,  and  z  frames,  DCR  is  modeled  as  a  linear 
function  of  color  (r—i)  for  all  stars  ([r— z]0  =  oo).  The  cor¬ 
rected  frame  coordinates  (x',  /)  are  then  transformed  to 
CMP  great  circle  coordinates  (/Amp?  *Amp)  using  an  affine 
transformation. 

Figure  2  displays  cubic  polynomial  fits  to  the  observed 
optical  distortions  as  a  function  of  CCD  column  for  each 
CCD.  Distortions  are  typically  less  than  0.2  pixels  and  are 
stable  from  month  to  month. 

The  use  of  a  set  of  transformation  equations  per  frame, 
rather  than  directly  calibrating  each  object,  facilitates  com- 
partmentalization  in  data  processing  as  well  as  later  recali- 
brations.  However,  it  can  introduce  systematic  errors  if  the 
equations  do  not  adequately  map  the  true  behavior.  After 
fitting  the  optical  distortion  terms,  there  remain  systematics 
as  a  function  of  column,  although  typically  of  the  order  of  5 
mas  or  less.  It  is  more  difficult  to  estimate  the  errors  intro¬ 
duced  by  the  approximation  that  the  astrometry  changes 
linearly  with  time  (that  is,  with  row)  over  a  single  frame  on 
the  photometric  CCDs;  however,  they  seem  to  be  less  than 
of  the  order  of  10  mas.  (This  is  not  a  good  approximation 
for  the  astrometric  CCDs,  for  which  the  effects  of  seeing  are 
more  severe  because  of  the  shorter  integration  times.) 

The  equations  for  fitting  DCR  were  determined  by  con¬ 
volving  spectrophotometry  of  stars  from  the  Gunn  & 
Stryker  (1983)  spectral  atlas  with  the  SDSS  filter  curves  and 


a  model  atmosphere.  For  the  r,  z,  and  z  filters,  the  linear  fits 
of  DCR  versus  color  at  a  given  air  mass  are  excellent  over 
all  spectral  types,  with  rms  errors  of  less  than  5  mas  at  a 
zenith  distance  of  60°  (the  extreme  zenith  distance  at  which 
observations  are  obtained;  as  zenith  distance  is  varied,  the 
DCR  errors  scale  with  refraction).  In  u  and  g  the  linear  fits 
are  poorer  but  adequate  if  confined  to  stars  bluer  than 
r—i  ~  1.5,  with  rms  errors  of  38  and  13  mas,  respectively,  at 
a  zenith  distance  of  60°.  For  redder  stars  the  DCR  terms  in 
u  and  g  are  poorly  behaved,  and  thus  a  constant  value  is 
employed.  The  accuracy  of  the  astrometry,  including  the 
systematics  introduced  by  the  form  of  the  transformation 
equations,  is  discussed  in  §  6. 

5.4.  r  Calibrations 

The  technique  used  to  calibrate  the  r  CCDs  depends  on 
the  primary  standard  star  catalog  used.  The  preferred  cata¬ 
log  is  UCAC,2  an  (eventually)  all-sky  astrometric  catalog 
with  a  precision  of  70  mas  at  its  catalog  limit  of  R  ~  16  and 
systematic  errors  of  less  than  30  mas.  Currently,  the  entire 
southern  hemisphere  has  been  observed,  with  partial  cover¬ 
age  in  the  north  (the  entire  sky  should  be  finished  by  the  end 
of  2003);  one-third  of  the  SDSS  area  is  covered  at  present. 
There  are  2-3  mag  of  overlap  between  the  faint  end  of 
UCAC  and  bright  unsaturated  stars  on  the  r  CCDs 
(depending  on  the  CCD  and  on  seeing),  allowing  the  direct 
calibration  of  the  r  CCDs  against  UCAC,  bypassing  the 
astrometric  CCD  to  photometric  CCD  bootstrap. 

If  a  scan  is  not  covered  by  the  current  UCAC  catalog, 
then  it  is  reduced  against  Tycho-2,  an  all-sky  astrometric 
catalog  with  a  median  precision  of  70  mas  at  its  catalog  limit 
of  Vj  ~  11.5  and  systematic  errors  of  less  than  1  mas.  All 
Tycho-2  stars  are  saturated  on  the  r  CCDs;  however,  there 
are  about  3.5  mag  of  overlap  between  bright  unsaturated 
stars  on  the  astrometric  CCDs  and  the  faint  end  of  Tycho-2 
(8  <  r  <  11.5)  and  about  3  mag  of  overlap  between  bright 
unsaturated  stars  on  the  r  CCDs  and  faint  stars  on  the  astro- 
metric  CCDs  (14<r<  17). For  reductions  against  Tycho- 
2,  the  star  centroids  on  the  astrometric  CCDs  are  trans¬ 
formed  to  the  pixel  coordinates  of  the  r  CCD  in  the  same 
column  using  the  unsaturated  stars  detected  on  both  the 
astrometric  and  r  CCDs.  Each  pair  of  astrometric -r  CCDs 
is  processed  separately.  Stars  on  each  astrometric  CCD  are 
matched  to  stars  on  their  associated  r  CCD.  There  are  typi¬ 
cally  10-20  matched  pairs  per  frame.  For  each  frame,  the 
nearest  (in  ju)  30  matched  pairs,  including  all  matched  pairs 
on  that  frame,  are  used  to  derive  an  affine  transformation 
from  the  astrometric  CCD  pixel  coordinate  system  to  the  r 
CCD  pixel  coordinate  system,  minimizing  the  rms  row  and 
column  differences  between  the  astrometric  and  r  stars.  If 
this  is  the  first  fitting  iteration,  outliers  are  rejected.  Smooth¬ 
ing  splines,  using  three  piecewise  polynomials  per  frame,  are 
then  fitted  to  the  row  and  column  residuals  as  a  function  of 
row  (time).  Outliers  are  rejected,  and  the  fitting  process  is 
iterated  two  more  times.  Finally,  the  remaining  residuals  for 
all  matched  pairs  are  fitted,  using  a  least-squares  fit  to  a 
cubic  polynomial,  as  a  function  of  column  to  remove  the 
optical  distortions.  The  resulting  solution  is  used  to  trans¬ 
form  the  stars  detected  on  the  astrometric  CCD  to  the  pixel 
coordinate  system  of  the  r  CCD.  If  a  star  was  detected  on 


2  See  http://ad.usno.navy.mil/ucac. 
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Fig.  2. — Cubic  polynomial  fits  to  the  observed  optical  distortions  for  TDI  scanning,  in  pixels,  as  a  function  of  CCD  column.  Plots  are  displayed  for  each 
CCD  (labeled  by  filter)  in  camera  columns  1-6  ( left  to  right )  of  the  imaging  camera.  The  distortions  along  the  rows  are  plotted  as  solid  lines  and  along  the 
columns  as  dotted  lines. 


both  the  leading  and  trailing  astrometric  CCDs  within  a  col¬ 
umn,  then  its  transformed  r  coordinates  are  averaged  and 
the  mean  position  is  used.  These  transformed  stars  will  then 
be  matched  to  the  Tycho-2  stars  to  yield  a  transformation 
from  r  pixel  coordinates  to  great  circle  CMP  coordinates. 

Figure  3  shows  a  typical  smoothing  spline  fit  to  the  row 
and  column  residuals  for  an  astrometric  CCD  to  r  CCD 
match.  The  difficulty  in  doing  astrometry  with  11s  exposures 
in  drift-scan  mode  is  evident.  Peak-to-peak  excursions  are 
100-300  mas,  with  timescales  of  the  order  of  half  a  frame  (18 
s).  These  most  likely  are  due  to  anomalous  refraction  (see 
§  6.4).  The  scan  used  for  the  plot  has  a  higher  number  of 
pairs  per  frame  than  most  scans,  with  half  that  number  more 
typical.  While  no  all- sky  astrometric  catalog  of  sufficient 
accuracy  is  dense  enough  to  follow  the  atmospheric  excur¬ 
sions  on  these  timescales,  there  are  enough  overlap  stars 
between  the  astrometric  and  r  CCDs  to  allow  most  of  the 
systematics  to  be  removed.  There  is  a  tendency  for  the  peaks 


of  the  excursions  to  be  underfitted.  Interpolating  splines 
were  experimented  with  in  place  of  the  smoothing  splines, 
but  the  final  astrometry  was  not  improved.  The  rms  residual 
in  either  row  or  column  for  the  transformation  from  an 
astrometric  CCD  to  an  r  CCD  is  typically  30-50  mas. 

The  reduction  now  proceeds  similarly  whether  reducing 
against  UCAC  or  Tycho-2.  The  following  focal  plane  model 
is  used,  relating  frame  row  (x)  and  column  (y)  to  PCP  coor¬ 
dinates  (the  small-angle  approximation  is  assumed  for  all 
rotations): 

x'  =  x  +  go  +  giy  +  g2y2  +  g2y3  ,  (11) 

/  =y  +  h0  +  hly  +  h2y2 +h2y3  ,  (12) 

/^*,PCP  =  /ffi,  CMP  +  T^x 

+  s  [-Tccd  —  Xb  +  (  Tccd  -  Yb)0 cam 
+  0.024/ccd (/  -  1024)(6>cam  +  #ccd)]  ,  (13) 
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Fig.  3. — Typical  smoothing  spline  fits  for  the  transformation  from  an 
astrometric  CCD  to  an  r  CCD.  The  20  frames  displayed  are  approximately 
12  minutes  of  data. 


^*,PCP  =  Vb,  CMP  +  Tvx' 

+  s[  —  (XcCD  -  Xb)0 cam  +  Yccd  -  YB 
+  0.024/ccd (y  —  1024)]  ,  (14) 

where  x'  and  /  are  frame  row  and  column  corrected  for 
optical  distortions.  The  telescope  tracking  (along  the  CMP 
great  circle)  is  modeled  by  specifying  the  boresight  pointing 
when  row  0  of  the  frame  was  read  (jib, cmp  and  +8, cmp?  in 
arcseconds)  and  the  tracking  rate  parallel  (7^,  in  arcsec 
pixel-1)  and  perpendicular  {Tv)  to  the  great  circle.  A  per¬ 
fectly  tracked  scan  would  have  vB^ cmp  =  0,  Tv  =  0,  and 
equal  to  the  sidereal  rate.  The  location  of  the  telescope  bore- 
sight  in  the  focal  plane  is  specified  by  the  constant  parame¬ 
ters  XB  (always  set  to  0)  and  YB  (—22.74  mm  for  north 
strips,  +22.74  mm  for  south  strips),  s  is  the  telescope  scale 
(in  arcsec  mm-1),  Oqam  is  the  rotation  of  the  camera  from 
its  correct  value  (in  radians),  Xqcd  and  Tccd  are  the  loca¬ 
tion  of  the  center  of  the  CCD  (column  1024)  in  the  focal 
plane  (in  millimeters),  Oqcd  is  the  rotation  of  the  CCD  with 
respect  to  the  focal  plane  coordinate  system, /ccd  is  the  rela¬ 
tive  scale  of  the  CCD  (with  a  nominal  value  of  1),  and  the 
size  of  the  pixels  is  0.024  mm.  The  quantities  PCp  and 
^*,pcp  are  the  resulting  PCP  coordinates  of  the  star  (in 
arcseconds). 

The  camera  geometry  (the  telescope  scale  and  the  CCD 
focal  plane  locations,  rotations,  relative  scales,  and  optical 
distortion  terms)  is  calibrated  once  a  month  using  a  scan 
through  a  UCAC  field.  (Before  UCAC  became  available, 
calibration  regions  measured  using  the  Flagstaff  Astrome¬ 
tric  Scanning  Transit  Telescope  [Stone  1997;  Stone,  Pier,  & 
Monet  1999]  were  used.)  The  most  recent  set  of  these  nomi¬ 
nal  values  is  used  as  an  initial  (and  in  some  cases  final)  guess 
at  those  values  for  the  current  reduction.  The  CCD  terms 
are  stable  from  month  to  month. 

The  focal  plane  model,  initialized  using  the  TCC  tele¬ 
scope  tracking  information  and  a  recent  calibration  scan,  is 


first  used  to  calculate  a  priori  CMP  coordinates  for  each 
observed  star  (stars  detected  on  the  r  CCDs  for  reductions 
against  UCAC  and  stars  detected  on  the  astrometric  CCDs 
and  converted  to  the  r  pixel  coordinate  systems  for  reduc¬ 
tions  against  Tycho-2).  These  a  priori  coordinates  are  used 
to  match  the  observed  stars  with  the  catalog  stars  (catalog 
stars  whose  catalog  error  in  either  right  ascension  or  decli¬ 
nation  exceeds  60  mas  are  not  used,  rejecting  about  25%  of 
the  stars).  For  each  matched  pair,  the  position  of  the  catalog 
star  is  converted  from  CMP  to  PCP  using  the  TT  at  mid¬ 
exposure  for  the  observed  star  and  the  color  of  the  star. 
UCAC  does  not  contain  color  information,  so  the  observed 
color  (r—i)  of  the  star  is  used  for  reductions  against  UCAC. 
Tycho-2  stars  saturate  the  photometric  CCDs  and  there¬ 
fore  do  not  have  SDSS  colors;  thus,  the  Tycho-2  catalog 
colors  (Bj—Vj)  are  used,  converted  to  r—i  using  the 
transformation 

r  —  i  —  —  0.186  +  0.681(£t  -  Ur) 

-  0.534(Rr  -  VT)2  +  0.319(_BT  -  Ft )3  •  (15) 

A  solution  to  the  focal  plane  model  is  now  derived  for 
each  frame  in  order  to  minimize  the  rms  differences  in  /q>cp 
and  v pep  between  the  observed  and  catalog  stars.  First,  a 
separate  linear  least-squares  fit  is  performed  for  each  set  of 
six  frames  observed  at  the  same  time  by  the  six  r  CCDs.  The 
CCD  terms  (focal  plane  locations,  scales,  rotations,  and 
optical  distortions)  are  held  constant  at  their  nominal 
values,  and  only  the  six  terms  specifying  the  boresight  track¬ 
ing  (jiBt cmp,  cmp?  Tp  and  Tv,  telescope  scale  [s],  and  cam¬ 
era  rotation  [#cam])  are  fitted.  For  a  given  set  of  six  frames, 
all  matched  pairs  on  that  set  of  frames  are  used.  Additional 
matched  pairs  from  the  nearest  adjacent  frames  are  used  if 
required  to  yield  a  minimum  number  of  matched  pairs  (30 
for  reductions  against  UCAC,  15  for  reductions  against 
Tycho-2).  If  an  adjacent  set  of  frames  is  used,  then  all 
matched  pairs  on  that  set  of  frames  are  used.  For  reductions 
against  UCAC,  there  are  typically  3-10  catalog  stars  per 
frame,  so  it  is  usually  not  necessary  to  use  the  adjacent 
frames.  For  reductions  against  Tycho-2,  there  is  typically 
one  catalog  star  every  two  frames;  thus,  it  is  generally  neces¬ 
sary  to  use  the  adjacent  four  or  more  sets  of  frames  to 
achieve  the  minimum  number  of  pairs. 

If  this  is  the  first  iteration,  outliers  are  rejected.  There  are 
systematics  remaining  in  the  residuals  after  the  least-squares 
fit,  primarily  because  of  anomalous  refraction,  varying  with 
typical  timescales  of  5-20  frames  (see  §  6.4).  For  reductions 
against  UCAC,  smoothing  splines  are  fitted  to  the  //  and  v 
residuals  for  each  r  CCD  separately,  with  breakpoints  sepa¬ 
rated  by  five  frames.  This  is  not  possible  for  reductions 
against  Tycho-2,  because  of  the  paucity  of  catalog  stars. 
Rather,  the  average  (over  the  entire  scan)  residual  in  n  and  v 
is  removed  for  each  r  CCD  separately.  In  either  case,  this 
procedure  effectively  refits  the  CCD  location  in  the  focal 
plane  (XCCD  and  FCCp).  For  both  the  UCAC  and  Tycho-2 
cases,  outliers  are  rejected,  and  the  process  of  the  least- 
squares  fit,  spline  fits,  and  outlier  rejection  is  iterated  two 
more  times. 

Figure  4  shows  typical  smoothing  spline  fits  to  the  fi  and  v 
residuals  for  all  six  r  CCDs  for  part  of  a  reduction  against 
UCAC  (separate  fits  for  each  CCD).  Systematics  in  com¬ 
mon  to  all  six  CCDs  (such  as  due  to  telescope  jitter  or  any 
portion  of  anomalous  refraction  that  is  correlated  across 
the  focal  plane)  have  already  been  removed  by  the  running 
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means.  The  remaining  systematics  have  typical  amplitudes 
of  20  mas,  with  timescales  of  a  few  frames  or  more.  There 
are  occasional  features  with  timescales  closer  to  a  frame, 
which  are  more  prevalent  in  worse  seeing.  The  spline  fits 
remove  most  of  the  systematics  but  not  all.  More  closely 
placed  breakpoints  have  been  experimented  with,  but  this 
can  lead  to  fitting  features  that  are  not  real.  The  chosen 
breakpoint  separation  of  five  frames  is  a  good  compromise 
between  fitting  most  of  the  real  features  without  overfitting 
the  data  for  typical  seeing  conditions.  These  systematics 
cannot  be  removed  for  reductions  against  Tycho-2,  as  there 
are  too  few  catalog  stars  detected  per  CCD  to  map  the  sys¬ 
tematics.  This  adds  roughly  10-20  mas  in  quadrature  to  the 
rms  errors  for  reductions  against  Tycho-2.  Similarly,  the 
running  mean  typically  smooths  over  single  frames  for 
reductions  against  UCAC  and  five  frames  for  reductions 
against  Tycho-2.  Thus,  any  correlated  systematics  with 
timescales  of  the  order  of  a  frame  cannot  be  removed  from 
reductions  against  Tycho-2. 

For  reductions  against  UCAC,  the  optical  distortion 
terms  (gb  hf)  are  then  refitted  by  fitting  the  remaining  resid¬ 
uals  to  a  cubic  polynomial  as  a  function  of  column.  One  set 
of  terms  is  calculated  per  CCD  for  the  entire  scan,  not  per 


frame.  This  effectively  refits  the  CCD  rotations  and  relative 
scales  ($ccd  and /ccd  in  the  focal  plane  model).  There  are 
not  enough  matched  pairs  for  a  reduction  against  Tycho-2, 
particularly  for  short  scans,  to  refit  these  terms,  so  the  nomi¬ 
nal  values  from  the  most  recent  calibration  scan  are  used. 
While  these  terms  are  very  stable  over  time,  occasionally 
there  is  evidence  of  residual  CCD  rotation  or  scale  changes. 

The  solution  is  a  transformation  (nonlinear  in  row)  from 
r  pixel  coordinates  to  PCP  coordinates.  This  needs  to  be 
converted  to  the  calibration  equations  of  §  5.3,  which  con¬ 
vert  from  pixel  coordinates  to  CMP  coordinates,  one  set  of 
equations  per  frame.  The  optical  distortion  terms  of  the  sol¬ 
ution  are  copied  directly  to  the  calibration  equations.  For 
each  frame,  PCP  coordinates  are  calculated  from  the  solu¬ 
tion  at  the  midpoint  of  each  of  the  four  sides  of  the  frame. 
These  are  converted  to  CMP  using  the  TT  at  midexposure 
and  assuming  a  star  of  zero  color.  These  four  test  points 
then  determine  the  affine  transformation  portion  of  the  cali¬ 
bration  equations.  The  CMP  is  similarly  calculated  for  11 
test  points,  all  located  in  the  middle  of  the  frame  but  span¬ 
ning  a  range  in  color.  The  /^cmp  and  z'cmp  offsets  from  a  star 
of  zero  color  are  then  fitted  as  linear  functions  of  color, 
yielding  the  DCR  terms  for  each  frame. 
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Fig.  4. — Typical  smoothing  spline  fits  for  the  transformation  from  r  CCDs  to  UCAC 
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5.5.  u,  g,  i,  andz  Calibrations 

The  r  calibration  equations  are  applied  to  all  the  stars 
detected  on  the  r  CCDs,  generating  CMP  positions  for  these 
stars.  This  creates  a  catalog  of  secondary  astrometric  stan¬ 
dards  that  is  used  to  calibrate  the  remaining  photometric 
CCDs.  Each  of  the  u ,  g ,  z,  and  z  CCDs  is  calibrated  sepa¬ 
rately.  The  stars  detected  on  a  given  CCD  are  matched 
against  the  secondary  catalog  from  the  r  CCD  in  the  same 
camera  column.  There  are  typically  10-40  matched  pairs 
per  frame  (5-20  on  the  u  CCDs).  For  each  matched  pair,  the 
position  of  the  catalog  star  is  converted  from  CMP  coordi¬ 
nates  to  PCP  coordinates  using  the  TT  at  midexposure  and 
the  color  (u—g  for  u,  g—r  for  g,  r—i  for  i  and  z)  of  the 
observed  star.  For  each  frame,  all  matched  pairs  on  that 
frame  are  used  to  derive  an  affine  transformation  from  pixel 
coordinates  to  PCP  coordinates,  minimizing  the  rms  differ¬ 
ences  in  //pep  and  z/PCP  between  the  observed  and  catalog 
stars.  If  there  are  fewer  than  15  matched  pairs  on  the  frame, 
then  the  nearest  matched  pairs  on  adjacent  frames  are  used 
to  guarantee  that  a  minimum  of  1 5  pairs  is  used. 

Outliers  are  rejected,  and  the  process  of  fitting  and  outlier 
rejection  is  iterated  two  more  times.  The  remaining  residuals 
for  all  matched  pairs  are  fitted  using  a  least-squares  fit  to  a 
cubic  polynomial  as  a  function  of  column  to  remove  the 
optical  distortions  (that  is,  one  polynomial  is  fitted  per  CCD 
for  the  entire  scan,  not  each  frame).  Chromatic  aberration 
in  the  optics  can  produce  systematics  with  color  in  the 
remaining  residuals.  This  is  significant  only  in  g.  Thus,  for  g 
only,  an  empirically  determined  linear  term  with  color 
(based  on  a  fit  to  a  single  scan,  but  which  has  proven  to  be 
stable  from  scan  to  scan)  is  removed  from  the  remaining 
residuals.  The  resulting  solution  transforms  pixel  coordi¬ 
nates  to  PCP  coordinates.  It  is  converted  to  the  calibration 
equations  of  §  5.3,  which  convert  from  pixel  coordinates  to 
CMP  coordinates,  using  the  same  procedure  used  for  the  r 
CCDs. 


6.  RESULTS 

This  section  quantifies  the  accuracy  of  the  astrometric 
calibrations.  The  calibrations  are  based  on  stars  brighter 
than  20th  magnitude,  for  which  centroiding  errors  are 
negligible  (except  in  the  u  filter,  in  which  the  centroiding 
errors  start  to  become  significant  by  20th  magnitude  and 
contribute  to  the  rms  errors  quoted  here).  Thus,  the  accu¬ 
racies  quoted  here  are  for  the  calibrations  only  and 
reflect  the  astrometric  accuracy  expected  for  stars 
brighter  than  20th  magnitude.  Centroiding  errors  (which 
are  calculated  by  Frames)  must  also  be  considered  for 
fainter  stars  and  all  galaxies  and  are  discussed  at  the  end 
of  the  section.  The  results  presented  here  are  based  on 
the  entire  SDSS  data  set  to  be  made  public  in  early  2003 
as  SDSS  data  release  1  (DR1).  This  includes  33  scans 
covering  roughly  1200  deg2  of  sky  (not  all  unique) 
reduced  against  UCAC  and  32  scans  covering  roughly 
750  deg2  of  sky  reduced  against  Tycho-2. 

The  distribution  of  internal  residuals  within  a  single  scan, 
the  distribution  of  differences  in  position  for  matched  pairs 
between  two  runs,  and  the  distribution  of  differences  in  posi¬ 
tion  for  the  SDSS  matched  against  the  external  catalogs 
used  in  this  analysis  are  all  well  characterized  by  Gaussians 
and  therefore  by  rms  residuals  and  differences.  Thus,  all 
errors  are  quoted  as  rms  residuals  or  differences  in  /i  and  v. 


6.1.  r  Astrometry 

The  absolute  accuracy  of  the  r  astrometry  is  difficult  to 
gauge,  as  there  are  no  astrometric  catalogs  as  deep  and 
accurate  as  the  SDSS  itself.  The  primary  internal  measure  of 
the  precision  of  the  astrometry  is  the  distribution  of  resid¬ 
uals  (measured  position  of  the  SDSS  object  minus  the  cata¬ 
log  position)  for  matched  pairs  (after  outlier  rejection)  on 
the  r  CCDs.  Figure  5  shows  histograms  of  these  residuals 
for  all  matched  pairs  for  reductions  against  UCAC  and 
Tycho-2  separately.  The  rms  residuals  per  coordinate  are  45 
and  75  mas  for  reductions  against  UCAC  and  Tycho-2, 
respectively. 

The  precision,  of  course,  can  vary  from  scan  to  scan  and 
during  individual  scans.  For  each  scan  the  rms  residuals 
(around  0)  in  //  and  v  are  calculated  separately  for  each 
frame,  for  each  r  CCD  separately.  To  assure  good  statistics, 
a  smoothing  window  is  used,  using  enough  adjacent  frames 
to  give  a  minimum  of  100  matched  pairs  per  frame.  These 
rms  residuals  in  /a  and  v  for  each  frame  are  recorded  with 
the  calibration  equations  for  that  frame  and  serve  as  the  best 
estimate  of  the  systematic  errors  for  that  frame.  Figure  6 
shows  histograms  of  these  rms  residuals  for  all  frames  in 
DR1.  All  the  distributions  peak  near  the  rms  values  indi¬ 
cated  in  Figure  5.  For  reductions  against  UCAC,  almost  all 
frames  have  rms  residuals  of  less  than  60  mas.  The  distribu¬ 
tions  for  reductions  against  Tycho-2  are  broader,  with  most 
frames  having  rms  residuals  of  less  than  100  mas.  The 
dependence  on  seeing  is  shown  in  Figure  7,  which  plots  the 
median  rms  residual  per  frame  binned  by  seeing  (as  mea¬ 
sured  by  PSP).  The  reductions  against  Tycho-2  are  consid¬ 
erably  more  sensitive  to  seeing  than  the  reductions  against 
UCAC,  demonstrating  that  the  longer  integration  times  on 
the  r  chips  (relative  to  the  astrometric  chips),  coupled  with 
the  greater  star  density  of  the  UCAC  catalog,  are  sufficient 
to  remove  most  of  the  atmospheric  effects. 

Some  parts  of  the  sky  have  been  scanned  multiple  times. 
Stars  in  nine  pairs  of  overlapping  scans  (including  some 


residual  (") 

Fig.  5. — Histograms  of  the  residuals  (SDSS  position  minus  catalog 
position)  for  scans  reduced  against  (a)  UCAC  and  (b)  Tycho-2.  The  solid 
and  dotted  histograms  are  for  the  ju  and  v  residuals,  respectively. 
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Fig.  6. — Histograms  of  the  rms  residuals  (SDSS  position  minus  catalog 
position)  for  each  frame  in  scans  reduced  against  (a)  UCAC  and  (b)  Tycho- 
2.  The  rms  residuals  for  each  frame  are  calculated  using  a  minimum  of  the 
nearest  100  matched  pairs.  The  solid  and  dotted  histograms  are  for  the  fi 
and  v  rms  residuals,  respectively. 

non-DR  1  scans),  reduced  against  both  UCAC  and  Tycho- 
2,  have  been  matched.  The  distribution  of  position  differen¬ 
ces  for  all  matched  pairs  is  shown  in  Figure  8.  The  rms  dif¬ 
ferences  are  about  50  mas  in  both  /i  and  v  for  reductions 
against  UCAC  and  73  mas  in  /i  and  85  mas  in  v  for  reduc¬ 
tions  against  Tycho-2.  These  are  smaller  than  one  might 
expect  from  simply  adding  the  internal  errors  in  quadrature. 
Errors  in  UCAC  and  Tycho-2  catalog  positions,  which 
propagate  to  SDSS  positions  calibrated  using  those  cata- 
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Fig.  7. — Median  rms  residual  per  frame  binned  by  seeing  for  the  r  CCDs. 
Triangles  are  for  the  n  residuals,  and  squares  the  v  residuals.  Open  symbols 
are  for  reductions  against  UCAC,  and  filled  symbols  are  for  reductions 
against  Tycho-2. 


Fig.  8. — Histograms  of  the  position  differences  for  matched  pairs  on 
overlapping  scans.  Results  from  reductions  against  UCAC  are  shown  in 
(a),  those  for  reductions  against  Tycho-2  in  (b).  The  solid  and  dotted  histo¬ 
grams  are  for  the  n  and  v  differences,  respectively. 

logs,  are  correlated  among  the  repeat  scans  and  thus  are 
partially  removed  when  differences  in  the  SDSS  positions 
are  computed. 

Variation  of  astrometric  precision  within  and  between 
scans  can  be  tested  by  sorting  the  matched  pairs  by  ti  and 
binning  them  into  groups  of  100,  separately  for  each  r  CCD. 
The  rms  differences  in  /i  and  v  are  then  calculated  for  each 
set  of  100  matched  pairs.  Figure  9  shows  histograms  of  these 
rms  differences.  For  reductions  against  UCAC,  the  distribu¬ 
tions  peak  around  45  mas,  with  most  of  the  data  having  rms 
differences  of  less  than  70  mas.  For  reductions  against 
Tycho-2,  the  distributions  are  fairly  broad,  with  typical  rms 
differences  of  60  mas  and  most  of  the  data  having  rms  differ¬ 
ences  of  less  than  130  mas. 

An  external  measure  of  the  accuracy  of  the  r  astrometry 
for  data  reduced  against  Tycho -2  can  be  obtained  by  re¬ 
reducing  those  scans  that  were  reduced  against  UCAC 
against  Tycho-2  and  matching  the  resulting  star  positions 
against  UCAC.  Histograms  of  the  position  differences  for 
all  matched  pairs  are  shown  in  Figure  10.  The  rms  differen¬ 
ces  are  76  mas  in  //  and  79  mas  in  v,  consistent  with  the  inter¬ 
nal  rms  residuals  (Fig.  5).  There  are  mean  offsets  of  —  15  mas 
in  //  and  18  mas  in  v.  These  are  likely  due  to  uncorrected 
magnitude  terms  in  star  positions  on  the  SDSS  astrometric 
CCDs,  due  to  charge  transfer  efficiency  (CTE)  effects  (seen 
by  matching  the  leading  and  trailing  astrometric  CCDs  in  a 
given  camera  column),  although  residual  magnitude  terms 
in  the  preliminary  UCAC  catalog  used  may  also  contribute. 
The  astrometric  CCD  magnitude  terms  are  not  currently 
corrected,  since  (1)  the  necessary  test  data  are  lacking,  and 
(2)  UCAC  coverage  should  be  full  sky  within  a  year,  allow¬ 
ing  all  SDSS  data  to  be  recalibrated  against  UCAC,  bypass¬ 
ing  the  astrometric  CCDs.  Again,  variation  of  astrometric 
accuracy  within  and  between  scans  can  be  tested  by  sorting 
the  matched  pairs  by  fi  and  binning  them  into  groups  of 
100,  separately  for  each  r  CCD.  The  rms  differences  in  /i 
and  v  are  then  calculated  for  each  set  of  100  matched  pairs. 
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Fig.  9. — Histograms  of  the  rms  position  differences  for  sets  of  100 
matched  pairs  between  overlapping  scans.  For  each  pair  of  overlapping 
scans  and  for  each  r  CCD,  the  matched  pairs  are  sorted  by  /i  and  binned 
into  groups  of  100  matched  pairs,  and  the  rms  differences  in  /u  and  v  are  cal¬ 
culated  for  each  set  of  100  matched  pairs.  Results  from  reductions  against 
UCAC  are  shown  in  ( a ),  and  those  for  reductions  against  Tycho-2  in  ( b ). 
The  solid  and  dotted  histograms  are  for  the  /i  and  v  rms  differences, 
respectively. 

Histograms  of  these  rms  differences  are  plotted  in  Figure  1 1 . 
The  distributions  are  similar  to  the  histograms  for  the  rms 
internal  residuals  (Fig.  6),  peaking  at  70  mas,  although  with 
a  somewhat  larger  fraction  of  the  data  with  rms  differences 
exceeding  100  mas. 

The  external  catalog  that  most  closely  matches  the  SDSS 
in  terms  of  depth,  accuracy,  and  sky  coverage  is  the  Two 
Micron  All  Sky  Survey  (2MASS;  Skrutskie  et  al.  1997). 


SDSSTycho_a  -  UCAC  (") 


Fig.  10. — Histograms  of  the  position  differences  for  matched  pairs  for 
reductions  against  Tycho-2  matched  against  UCAC.  The  solid  and  dotted 
histograms  are  for  the  n  and  v  differences,  respectively. 


Fig.  11. — Histograms  of  the  rms  position  differences  for  sets  of  100 
matched  pairs,  binned  in  //,  for  reductions  against  Tycho-2  matched  against 
UCAC.  For  each  scan  and  for  each  r  CCD,  the  matched  pairs  are  sorted  by 
/i  and  binned  into  groups  of  100  matched  pairs,  and  the  rms  differences  in  /i 
and  v  are  calculated  for  each  set  of  100  matched  pairs.  The  solid  and  dotted 
histograms  are  for  the  fi  and  v  rms  differences,  respectively. 

Typical  position  errors  for  the  2MASS  second  incremental 
release  are  quoted  as  about  140  mas,  although  those  are 
probably  overstated,  and  a  truer  measure  of  the  rms  errors 
may  be  closer  to  1 10  mas  (Cutri  et  al.  2000).  Figure  12  shows 
histograms  of  the  differences  between  the  SDSS  and  2MASS 
positions  for  stars  in  nine  scans.  2MASS  sources  were 


Fig.  12. — Histograms  of  the  position  differences  for  matched  pairs 
between  the  SDSS  and  2MASS  for  nine  scans.  Results  from  reductions 
against  UCAC  are  shown  in  (< a ),  and  those  for  reductions  against  Tycho-2 
in  ( b ).  The  solid  and  dotted  histograms  are  for  the  //  and  v  differences, 
respectively. 
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limited  to  nonconfused,  nonextended  sources  detected  in  all 
three  bands  with  J  <  15.8,  H  <  15.1,  and  K  <  14.3  and 
with  the  major  axis  of  the  position  error  ellipse  less  than  175 
mas.  The  rms  differences  are  118  mas  in  /a  and  114  mas  in  v 
for  SDSS  data  reduced  against  UCAC  and  123  mas  in  /i  and 
122  mas  in  v  for  SDSS  data  reduced  against  Tycho-2.  These 
values  are  roughly  what  one  would  expect  by  combining  the 
rms  errors  in  the  two  catalogs  in  quadrature.  There  are 
mean  offsets  of  —25  mas  in  /a  and  —58  mas  in  v  for  SDSS 
reductions  against  UCAC  and  -35  mas  in  //  and  -33  mas  in 
v  for  SDSS  reductions  against  Tycho-2.  While  the  offsets  in 
//  could  be  explained  by  systematic  errors  in  the  SDSS 
reductions,  the  v  offsets  cannot.  Similar  offsets  in  v  are  seen 
when  matching  UCAC  against  2MASS  in  the  same  region 
of  sky.  A  known  bias  exists  in  the  2MASS  second  incremen¬ 
tal  release  reductions  whereby  the  positions  of  bright  stars 
(K<  9),  which  are  processed  through  the  Rl  path  of  the 
2MASS  pipeline,  can  be  offset  from  those  of  fainter  stars, 
which  are  processed  through  the  Rl  -  Rl  path  of  the 
2MASS  pipeline;  this  bias  will  be  addressed  in  the  next 
2MASS  data  release  (H.  L.  McCallon  2001,  private  commu¬ 
nication).  This  bias  is  likely  the  largest  contributor  to  the  v 
offsets  between  the  SDSS  and  2MASS  and  may  also  contrib¬ 


ute  to  the  n  offsets.  (An  earlier  comparison  between  SDSS 
and  2MASS  astrometry,  based  on  a  smaller  data  set,  is  pre¬ 
sented  by  Finlator  et  al.  2000.) 

It  should  be  emphasized  that  the  values  of  the  rms  resid¬ 
uals  and  differences  quoted  characterize  the  distribution  of 
systematic  errors.  The  analysis  is  based  on  stars  with  r  <  20, 
for  which  random  centroiding  errors  are  small  compared  to 
the  systematic  errors.  Figure  13  displays  the  differences  in  fi 
and  v  between  matched  stars  in  two  overlapping  scans 
reduced  against  UCAC  plotted  against  fi  (time).  The  /i  and 
v  differences  are  plotted  separately  for  each  r  CCD. 
Figure  14  shows  the  differences  for  matches  for  the  same 
portion  of  the  same  overlapping  scans  but  reduced  against 
Tycho-2.  The  differences  vary  systematically  with  time,  on 
timescales  of  the  order  of  1  to  a  few  minutes.  While  the  dif¬ 
ferences  are  mostly  uncorrelated  between  CCDs  for  the 
reductions  against  UCAC,  they  are  strongly  correlated  for 
the  reductions  against  Tycho-2.  There  are  simply  too  few 
Tycho-2  stars  to  adequately  follow  the  tracking  errors  and 
that  portion  of  anomalous  refraction  that  is  correlated 
across  the  camera  on  short  timescales.  These  plots  are  typi¬ 
cal  and  display  well  the  nature  of  the  dominant  systematics 
in  the  astrometry. 
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Fig.  13. — Differences  in  position  between  matched  stars  in  two  overlapping  scans  reduced  against  UCAC  plotted  against  //  (time).  The  fi  and  v  differences 
are  plotted  separately  for  each  r  CCD. 
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Fig.  14. — Same  as  Fig.  13,  but  with  data  reduced  against  Tycho-2 


SDSS  astrometry  is  subject  to  any  systematics  present  in 
the  primary  astrometric  reference  catalog.  At  the  current 
epoch,  systematics  in  Tycho-2  are  expected  to  be  less  than  a 
few  milliarcseconds,  and  as  such  do  not  degrade  the  overall 
accuracy  of  SDSS  astrometry.  Uncorrected  magnitude 
terms  in  the  star  positions  on  the  astrometric  CCDs,  used  to 
transfer  the  Tycho-2  stars  to  the  photometric  CCDs,  intro¬ 
duce  systematic  offsets  between  the  Tycho-2  frame  and 
SDSS  data  reduced  against  Tycho-2  of  the  order  of  20  mas. 
A  preliminary,  internal  version  of  UCAC  (r07)  is  currently 
being  used  for  SDSS  astrometric  reductions.  This  version  of 
the  catalog  contains  some  uncorrected  magnitude  terms, 
such  that  the  brighter  stars,  which  are  used  to  calibrate 
UCAC  against  Tycho-2,  may  be  systematically  offset  from 
the  fainter  stars,  which  are  used  to  calibrate  the  SDSS.  Such 
offsets  are  thought  to  be  of  the  order  of  10-20  mas,  mainly 
caused  by  residual  CTE  effects  in  the  CCD  used  for  the 
UCAC  observations.  In  addition,  there  will  be  zero-point 
errors  due  to  the  fact  that  a  given  UCAC  frame  can  only  be 
linked  with  a  certain  accuracy  to  the  Tycho-2  stars,  depend¬ 
ing  on  the  number  of  reference  stars  available  in  that  area  of 
the  sky.  This  will  cause  systematic  errors  of  up  to  30  mas  on 
the  0?5-l°  scale.  Since  UCAC  uses  only  a  single  bandpass, 
colors  for  the  stars  were  not  available  and  DCR  corrections 


could  not  be  applied.  However,  the  narrow  bandpass  used 
causes  only  5-10  mas  of  systematic  errors  as  a  function  of 
spectral  type  and  zenith  distance.  In  sum,  total  systematic 
errors  in  UCAC  are  thought  not  to  exceed  30  mas  (N. 
Zacharias  2001,  private  communication).  As  newer  versions 
of  UCAC  with  smaller  systematic  errors  become  available, 
the  SDSS  astrometry  will  be  recalibrated. 

6.2.  u,  g,  i,  and  z  Astrometry 

The  primary  internal  measure  of  the  accuracy  of  the  rela¬ 
tive  astrometry  between  filters  is  the  distribution  of  residuals 
(measured  positions  in  the  target  filter,  u ,  g,  z,  or  z,  minus 
the  measured  positions  in  the  r  filter).  Figure  1 5  shows  histo¬ 
grams  of  these  residuals  for  all  matched  pairs  (the  relative 
astrometry  between  filters  is  not  sensitive  to  the  choice  of 
primary  reference  catalog),  separately  for  each  filter.  The 
relative  astrometry  is  best  in  the  g  and  z  filters,  with  rms 
residuals  of  about  24  mas  per  coordinate.  The  z  filter  is  only 
slightly  worse,  with  rms  residuals  of  about  27  mas  per  coor¬ 
dinate.  The  centroids  are  poorer  in  the  u  filter  because  of  the 
lower  flux  for  the  stars  used  to  derive  the  transformations, 
leading  to  larger  rms  residuals  of  about  33  mas  per  coordi¬ 
nate.  Again,  the  variation  in  astrometric  accuracy  can  be 
examined  by  calculating  rms  values  for  each  frame,  for  each 
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Fig.  15. — Histograms  of  the  residuals  (SDSS  position  in  the  specified  filter  minus  SDSS  position  in  the  r  filter)  for  all  DR1  scans.  The  solid  and  dotted 
histograms  are  for  the  fi  and  v  residuals,  respectively. 


CCD  pair.  To  assure  good  statistics,  a  smoothing  window  is 
used,  using  enough  adjacent  frames  to  give  a  minimum  of 
100  matched  pairs  per  frame.  The  rms  residuals  in  n  and  v 
for  each  frame  are  recorded  with  the  calibration  equations 
for  that  frame  and  serve  as  the  best  estimate  of  the  system¬ 
atic  error  for  that  frame.  Figure  16  shows  histograms  of 
these  rms  residuals.  The  distributions  are  fairly  narrow  in  g, 
z,  and  z,  with  most  frames  having  rms  residuals  of  less  than 
35  mas.  The  distribution  in  u  shows  a  tail  out  to  60  mas.  This 
is  due  to  both  the  increased  centroiding  errors  in  the  stars 
used  to  calibrate  the  u  frames  and  an  occasional  lack  of 
enough  bright  stars  to  adequately  calibrate  the  frame.  The 
rms  residuals  are  sensitive  to  seeing,  as  shown  in  Figure  17, 
which  plots  the  median  rms  residual  per  frame  binned  by 
seeing. 

Figures  18-20  plot  the  //  and  v  residuals  against  CCD 
column,  magnitude,  and  color,  respectively,  for  all  five 
photometric  CCDs  in  column  3  of  the  camera  for  a  scan 
taken  under  exceptional  seeing  to  minimize  the  random 
errors  and  thus  make  any  remaining  systematic  errors 
more  obvious.  The  scatter  for  the  r  CCD  is  larger  than 
the  other  CCDs  because  r  has  been  matched  against 
UCAC  while  the  other  CCDs  were  matched  against  the  r 
detections.  Systematics  with  column  (uncorrected  by  the 
cubic  polynomials  used  to  remove  the  optical  distortions) 
are  almost  always  less  than  10  mas  (although  for  some  u 
CCDs  they  can  peak  at  15  mas)  and  typically  are  5  mas 
or  less.  Systematics  with  magnitude  are  always  less  than 
13  mas  over  the  magnitude  range  14  <  r  <  22  and  are 
typically  less  than  half  that.  Systematics  with  color  are 
evident  only  for  the  u  and  g  CCDs  and  are  typically  less 
than  10  mas  over  the  full  range  in  color.  All  of  these 


uncorrected  systematics  contribute  to  the  rms  residuals 
quoted  above. 

6.3.  Centroiding  Errors 

The  results  presented  above  apply  to  stars  brighter  than 
r  ~  20,  for  which  systematic  errors  dominate  centroiding 
errors  and  thus  characterize  the  errors  in  the  calibrations. 
Centroiding  errors  become  significant  for  stars  fainter  than 
r  ~  20  and  for  all  galaxies.  Figure  21  shows,  in  the  top  two 
panels,  the  position  differences  against  r  for  matched  stars 
between  two  overlapping  scans  reduced  against  UCAC. 
Each  scan  was  taken  under  fairly  constant  seeing  of  1  '/4  (the 
median  seeing  for  all  DR1  data  is  about  l'/4).  The  errors 
clearly  increase  for  stars  fainter  than  r  ~  20.  The  position 
errors  as  a  function  of  magnitude  for  a  run  taken  during  see¬ 
ing  of  1  '!4  can  be  estimated  by  binning  the  position  differen¬ 
ces  in  magnitude,  calculating  the  rms  differences  in  each 
magnitude  bin,  and  dividing  the  rms  by  \fl.  This  has  been 
done  in  the  bottom  panels  of  Figure  21.  The  rms  errors 
increase  from  about  40  mas  at  r  ~  20  to  about  100  mas  at 
the  survey  limit  of  r  ~  22.  The  errors  are  seeing  dependent, 
with  rms  errors  at  r  ~  22  ranging  from  about  70  to  140  mas 
in  seeing  of  approximately  1  "A  and  1  ''7,  respectively  (a  range 
of  seeing  that  characterizes  most  of  the  DR1  data). 

Figure  22  is  similar  to  Figure  21  but  plots  the  position  dif¬ 
ferences  and  estimated  errors  for  galaxies  rather  than  stars 
(for  the  same  scans).  Centroiding  errors  dominate  the  cali¬ 
bration  errors  even  at  bright  magnitudes,  with  the  rms 
errors  ranging  from  about  60  mas  at  r  ~  17  to  170  mas  at 
r  ~  22.  These  results  are  independent  of  seeing. 

The  accuracy  of  the  relative  astrometry  between  filters 
similarly  decreases  with  increasing  magnitude.  Figure  23 


Fig.  16. — Histograms  of  the  rms  residuals  (SDSS  position  in  the  specified  filter  minus  SDSS  position  in  the  r  filter)  for  each  frame  in  all  DR1  scans.  The  rms 
residuals  for  each  frame  are  calculated  using  a  minimum  of  the  nearest  100  matched  pairs.  The  solid  and  dotted  histograms  are  for  the  //  and  v  rms  residuals, 
respectively. 


Fig.  17. — Median  rms  residual  per  frame  binned  by  seeing  for  the  u,  g,  i,  and  z  CCDs.  Triangles  are  for  the  //  residuals,  and  squares  the  v  residuals. 
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Fig.  19. — Same  as  Fig.  1 8,  but  plotted  against  magnitude 
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Fig.  20. — Same  as  Fig.  1 8,  but  plotted  against  color  (r—i  for  r,  i,  and  z,  u—g  for  u ,  and  g—r  for  g) 
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Fig.  23. — Position  differences  between  the  i,  u,  z,  and  g  filters  and  the  r  filter,  plotted  against  magnitude  (i,  u,  z,  g)  for  a  typical  scan 
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Fig.  24. — Plots  of  nns  errors  for  the  relative  astrometry  between  the  i,  u,  z,  and  g  filters  and  the  r  filter,  plotted  against  magnitude  (;,  u,  z,  g)  for  a 
typical  scan. 


plots  the  position  differences  between  the  z,  u ,  z,  and  g  filters 
and  the  r  filter  against  magnitude  (z,  u,  z ,  g )  for  a  single  scan 
observed  in  l'/65  seeing.  The  rms  errors  as  a  function  of 
magnitude  can  be  calculated  for  this  scan  by  binning  the  dif¬ 
ferences  in  magnitude  and  calculating  the  rms  differences  in 
each  magnitude  bin.  These  errors  are  shown  in  Figure  24. 
The  rms  errors  increase  to  around  100-130  mas  by  the 
survey  limits. 

6.4.  Anomalous  Refraction 

The  primary  contributions  to  systematic  astrometric 
residuals  appear  to  come  from  one  or  both  of  two  causes: 
atmospheric  fluctuations  and/or  telescope  jitter.  These 
observed  systematics  display  considerable  coherence  across 
the  focal  plane;  that  is,  both  the  direction  and  amplitude  of 
the  residuals  from  CCD  to  CCD  show  a  high  degree  of  cor¬ 
relation.  These  effects  can  be  seen  in  Figures  13  and  14.  The 
focal  plane  subtends  2?3  across  the  sky,  and  the  physical 
dimension  is  on  the  order  of  50  cm  (i.e.,  several  times  typical 
values  for  r0). 

Thus,  one  is  led  to  believe  that  either  the  focal  plane  itself 
is  moving  (due,  e.g.,  to  telescope  jitter)  or  the  atmosphere  is 
coherently  moving  the  telescope  beam.  These  two  effects  are 
difficult  to  disentangle,  and  it  is  with  some  caution  that  we 
assign  these  effects  to  the  atmosphere  for  the  following 
reasons: 

Among  the  factors  one  might  imagine  could  cause  the 
telescope  to  be  the  culprit  are  mechanical  resonances,  servo 
motor  errors,  and/or  bearing  unevenness,  which  result  in 
tracking  errors,  wind-induced  motions,  thermal  oscilla¬ 
tions,  and  oscillations  in  the  mirror  support  system.  Time- 
scales  for  the  systematic  effects  vary  from  night  to  night 


(and,  occasionally,  from  hour  to  hour)  but  are  typically  on 
the  order  of  a  few  minutes  to  a  few  tens  of  minutes.  The  tele¬ 
scope  mechanical  resonances  are  on  the  order  of  a  few  hertz 
or  higher  and  would  average  out  over  an  integration  time. 
These  same  effects  have  been  seen  on  engineering  scans  dur¬ 
ing  which  the  telescope  was  clamped  on  the  meridian  and 
equator  with  all  telescope  drive  motors  turned  off.  During 
such  observations,  no  significant  motion  of  either  the  tele¬ 
scope  mechanicals  or  optics,  or  of  the  imaging  camera  and/ 
or  its  CCDs,  is  thought  to  have  contributed  to  the  effects 
during  these  scans.  Hence,  we  are  persuaded  that  neither 
tracking  errors  nor  mirror  motions  cause  these  effects.  Ther¬ 
mal  timescales  are  much  longer  (an  hour  or  more)  and  thus 
also  seem  to  be  an  unlikely  cause.  Wind-induced  motions 
remain  a  possibility,  although  we  see  such  effects  on  all  of 
our  scans,  no  matter  the  wind  speed  or  direction,  and  the 
telescope  is  enshrouded  in  a  mechanically  isolated  wind 
baffle  to  prevent  just  such  a  problem. 

Furthermore,  this  kind  of  behavior  has  been  seen  before 
in  drift-scanning  CCD  observations.  Dunham,  McDonald, 
&  Elliot  (1991)  observed  several  strips  multiple  times  with  a 
scanning  CCD.  A  comparison  of  the  residuals  “revealed 
the  presence  of  a  pervasive  low-frequency  motion  of  the  sky 
coordinate  system  relative  to  the  CCD.”  They  also  found 
that  the  variations  along  right  ascension  and  declination 
were  significantly  different.  Stone  et  al.  (1996)  reported  simi¬ 
lar  results  with  their  drift-scan  mode  CCD  observations 
made  on  an  extremely  stable  transit  telescope.  Both  of  these 
sets  of  observations  were  obtained  with  single  CCDs  (a 
couple  of  centimeters  in  size  at  most)  with  fields  of  view  of 
about  20'.  Stone  et  al.  (1996)  attribute  the  cause  of  these 
effects  to  “  anomalous  refraction,”  i.e.,  refraction  that  varies 
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from  the  smooth  analytical  models  of  refraction  that  are 
functions  only  of  zenith  distance. 

The  SDSS  observations  are,  we  believe,  the  first  astro¬ 
nomical  observations  to  demonstrate  that  these  effects  are 
due  to  phenomena  whose  scale  is  on  the  order  of  2°  or  more 
on  the  sky.  Webb  (1984)  and  Naito  &  Sugawa  (1984)  discuss 
atmospheric  boundary  layers  with  heights  of  typically  a  few 
hundred  meters  to  2  km  above  the  earth’s  surface  and  hori¬ 
zontal  wavelengths  (A)  of  1-10  or  more  km,  which  may  give 
rise  to  anomalous  refraction.  If  this  is  true,  then  we  would 
expect  to  see  characteristic  timescales  of  X/v  (where  v  is  the 
wind  speed).  For  typical  values  (1  km  <  A  <  10  km  and 
10  km  hr-1  <  v  <  40  km  hr  *),  we  expect  to  see  time  signa¬ 
tures  of  a  few  to  several  tens  of  minutes,  in  agreement  with 
the  SDSS  observations. 

7.  EARLY  DATA  RELEASE 

This  paper  describes  the  astrometric  calibration  proce¬ 
dure  based  on  Astrom,  Version  3.6,  and  Frames,  Version 
5.3.  On  2001  June  5,  the  SDSS  made  its  first  public  data 
release  of  primarily  commissioning  data,  referred  to  as  the 
Early  Data  Release  (EDR;  Stoughton  et  al.  2002).  The 
EDR  data  were  processed  with  earlier  versions  of  both 
Astrom  and  Frames,  and  so  the  astrometric  calibration  pro¬ 
cedure  differed  to  some  extent  from  the  procedure  described 
in  this  paper.  The  primary  difference  for  astrometry  con¬ 
cerns  the  centroids,  which  were  not  corrected  for  asym¬ 
metries  in  the  PSFs. 

More  significantly,  Frames  used  a  different  smoothing 
length,  as  well  as  different  magnitude-dependent  binning 
factors,  than  did  SSC  (earlier  versions  of  Astrom  used 
uncorrected  centroids  from  SSC  rather  than  the  currently 
corrected  centroids  from  PSP)  before  calculating  the  cen¬ 
troids.  This  led  to  systematic  offsets  between  the  SSC  and 
Frames  centroids  of  up  to  20  mas,  as  well  as  systematic  off¬ 
sets  between  Frames  centroids  for  bright  and  faint  objects 
of  the  order  of  20  mas  (the  break  occurs  around  r  ~  20). 
Since  the  calibrations  were  based  on  SSC  centroids,  but  the 
final  positions  are  based  on  Frames  centroids,  this  contrib¬ 
uted  an  additional  20  mas  of  systematic  error  in  the  EDR 
data  for  bright  (r<20)  objects,  and  perhaps  twice  that  for 
faint  objects.  Roughly  90%  of  the  data  in  the  EDR  have 
been  recalibrated  using  the  new  procedure,  and  these  new 
calibrated  positions  will  be  included  in  DR1 . 

8.  SUMMARY 

The  astrometric  calibration  of  the  SDSS  uses  two  reduc¬ 
tion  strategies.  In  those  areas  of  the  sky  with  UCAC  cover¬ 
age,  the  r  CCDs  are  calibrated  directly  against  UCAC.  In 
those  areas  of  the  sky  that  lack  UCAC  coverage,  bright  stars 
detected  on  the  astrometric  CCDs  are  transferred  to  the 
pixel  coordinate  systems  of  the  r  CCDs  using  stars  in  com¬ 
mon  to  both  CCDs,  and  the  bright  stars  are  then  used  to  cal¬ 
ibrate  the  r  CCDs  against  Tycho-2.  The  accuracy  of  the 
calibrations  using  UCAC  are  of  the  order  of  45  mas  rms, 
with  an  additional  systematic  error  of  up  to  30  mas  (due  pri¬ 
marily  to  systematic  errors  in  UCAC).  The  accuracy  of  cali¬ 
brations  using  Tycho -2  are  of  the  order  of  75  mas  with  an 
additional  systematic  error  of  the  order  of  20  mas  (due  to 
CTE  effects  in  the  astrometric  CCDs).  The  rms  errors  are 
dominated  by  Gaussian  distributions  of  systematic  errors 
that  vary  on  timescales  of  1  to  a  few  tens  of  minutes  because 


of  anomalous  refraction  and  by  random  errors  in  the  pri¬ 
mary  reference  catalogs.  The  accuracy  of  the  relative 
astrometry  of  the  u ,  g ,  i,  and  z  filters  versus  the  r  filter  is  of 
the  order  of  25  mas  rms  for  the  g,  i,  and  z  filters  and  35  mas 
rms  for  the  u  filter.  Systematic  errors  with  magnitude,  color, 
or  CCD  column  are  typically  less  than  10  mas. 

The  astrometric  performance  being  realized  by  SDSS  far 
exceeds  the  original  design  specifications  and  science 
requirements.  There  are  a  number  of  reasons  for  this: 

1 .  The  imaging  camera  has  proven  to  be  rigid  and  stable 
and  makes  no  significant  contribution  to  astrometric  errors. 
It  was  originally  anticipated  that  Dewar-to-Dewar  excur¬ 
sions  of  a  few  microns  (the  scale  is  ^60  jum  arcsec-1)  and  an 
undetermined  although  small  amount  of  flexure  might 
need  to  be  tolerated.  The  camera’s  performance  has  been 
exceptional. 

2.  The  telescope  tracking  has  proven  to  be  extremely 
accurate  and  reliable.  Considerable  engineering  and  design 
effort  went  into  optimizing  the  SDSS  2.5  m  telescope  per¬ 
formance,  with  obvious  benefits  for  astrometry  and  imaging 
quality. 

3.  The  availability  of  new,  deeper,  and  better  astrometric 
catalogs  (e.g.,  Tycho-2  and  UCAC)  has  made  a  significant 
impact  on  improving  the  astrometry. 

UCAC  provides  a  valuable  densification  of  the  Hipparcos 
Reference  Frame  (HRF),  the  best  current  optical  realization 
of  the  International  Celestial  Reference  Frame  (ICRF), 
extending  it  to  R  ~  16.  Ultimately,  the  SDSS  will  provide, 
for  about  a  fourth  of  the  sky,  a  further  densification  of  the 
HRF,  providing  astrometry  relative  to  the  HRF  (via  its  cali¬ 
bration  against  UCAC)  accurate  to  45  mas  rms  plus  30  mas 
systematic  to  20,  and  approximately  100  mas  rms  at 
r  ~  22. 

The  accuracy  of  the  astrometry  relative  to  the  HRF 
allows  for  precise  matching  with  other  catalogs  that  have 
also  been  calibrated  within  the  HRF.  Various  proper- 
motion  studies  are  currently  underway  based  on  matches  to 
USNO-B  (Monet  et  al.  2003).  Matches  between  the  SDSS 
and  2MASS  have  been  used  to  study  the  optical  and  infra¬ 
red  photometric  properties  of  stars  (Finlator  et  al.  2000) 
and  to  search  for  z  >  5.8  quasars  (Fan  et  al.  2001).  Matches 
between  the  SDSS  and  the  Faint  Images  of  the  Radio  Sky  at 
Twenty  cm  (FIRST;  Becker,  White,  &  Helfand  1995)  survey 
are  being  used  to  study  the  optical  and  radio  properties  of 
extragalactic  sources  (Ivezic  et  al.  2002).  The  ease  with 
which  such  studies  can  be  made  is  directly  attributable  to 
the  impact  of  the  watershed  Hipparcos  space  astrometry 
mission  (Perryman  et  al.  1997)  on  global  astrometry.  With 
the  accuracy  of  the  relative  astrometry  between  filters,  aster¬ 
oids  are  easily  detected  based  on  their  relative  motions 
between  filters  within  a  single  scan  (Ivezic  et  al.  2001),  and 
the  same  technique  offers  promise  for  the  detection  of 
Kuiper  Belt  objects. 

Funding  for  the  creation  and  distribution  of  the  SDSS3 
archive  has  been  provided  by  the  Alfred  P.  Sloan  Founda¬ 
tion,  the  Participating  Institutions,  the  National  Aero¬ 
nautics  and  Space  Administration,  the  National  Science 
Foundation,  the  US  Department  of  Energy,  the  Japanese 
Monbukagakusho,  and  the  Max  Planck  Society. 


3  The  SDSS  Web  site  is  http://www.sdss.org. 
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The  SDSS  is  managed  by  the  Astrophysical  Research 
Consortium  (ARC)  for  the  Participating  Institutions.  The 
Participating  Institutions  are  the  University  of  Chicago, 
Fermilab,  the  Institute  for  Advanced  Study,  the  Japan 
Participation  Group,  The  Johns  Hopkins  University,  Los 
Alamos  National  Laboratory,  the  Max-Planck-Institute  for 
Astronomy,  the  Max-Planck-Institute  for  Astrophysics, 
New  Mexico  State  University,  the  University  of  Pittsburgh, 
Princeton  University,  the  US  Naval  Observatory,  and  the 
University  of  Washington. 

Fermilab  is  operated  by  Universities  Research  Associa¬ 
tion,  Inc.,  under  contract  DE-AC02-76CH03000  with  the 
US  Department  of  Energy.  R.  H.  L.  and  Z.  I.  would  like  to 
thank  Princeton  University  for  generously  providing 
research  funds  and  Jill  Knapp  for  generously  providing 
lunch  and  sympathy.  We  thank  the  entire  UCAC  team  for 
making  available  to  the  SDSS  prerelease  versions  of  UCAC, 


Mario  Nonino  for  work  on  an  early  version  of  the  pipeline, 
Walter  Siegmund  for  useful  discussions,  and  N.  Zacharias, 
D.  York,  and  M.  Strauss  for  critical  readings  of  the  manu¬ 
script.  Astrom  makes  extensive  use  of  the  SLALIB  library 
(Wallace  1994)  from  the  Starlink  Project,  which  is  run  by 
the  Council  for  the  Central  Laboratory  of  the  Research 
Councils  on  behalf  of  the  Particle  Physics  and  Astronomy 
Research  Council.  This  publication  makes  use  of  data  prod¬ 
ucts  from  2MASS,  which  is  a  joint  project  of  the  University 
of  Massachusetts  and  the  Infrared  Processing  and  Analysis 
Center,  California  Institute  of  Technology,  funded  by  the 
National  Aeronautics  and  Space  Administration  and  the 
National  Science  Foundation.  This  research  has  made  use 
of  the  NASA/IPAC  Infrared  Science  Archive,  which  is 
operated  by  the  Jet  Propulsion  Laboratory,  California  Insti¬ 
tute  of  Technology,  under  contract  with  the  National 
Aeronautics  and  Space  Administration. 
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